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AERONAUTIC  SYMBOLS 

I.  FUNDAMENTAL  AND  DERIVED  UNITS 


Symbol 

Metric 

English 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia¬ 

tion 

l 

t 

F 

meter _ _ _ 

m 

s 

kg 

foot  (or  mile) _ 

ft  (or  mi) 
seo  (or  hr) 
lb 

Force . 

weight  of  1  kilogram . 

weight  of  1  pound _ 

P 

V 

horsepower  (metric) _ 

horsepower _ 

bp 

mph 

fps 

Speed..  _ 

/kilometers  per  hour _ 

(meters  per  second _ 

kph 

mp3 

miles  per  hour _ 

feet  per  second _ 

w 

a 

m 

1 

S 

sm 

a 

b 

c 

A 

V 

2 
L 
D 
Dt 
D, 
D, 
a 


2.  GENERAL  SYMBOLS 


Weight =m# 

Standard  acceleration  of  gravity=9.80665  m/s1 
or  32.1740  ft/sec3 
W 

Mass=y 

Moment  of  inertia=mt*.  (Indicate  axis  of 
radius  of  gyration  k  by  proper  subscript.) 
Coefficient  of  viscosity 


,  Kinematic  viscosity 

p  Density  (mass  per  unit  volume) 

Standard  density  of  dry  air,  0.12497  kg-m-i-s*  at  15°  C 
and  760  mm;  or  0.002378  lb-ft-4  sec3 
Specific  weight  of  “standard”  air,  1.2255  kg/m*  or 
0.07651  lb/cu  ft' 


X  AERODYNAMIC  SYMBOLS 


Area 

Area  of  wing 
Gap 
Span 
Chord 

Aspect  ratio,  g- 
True  air  speed 
Dynamic  pressure,  |pV* 


lift,  absolute  coefficient  °L=fs 
Drag,  absolute  coefficient  CrD=^ 
Profile  drag,  absolute  coefficient 


Induced  drag,  absolute  coefficient  Cd<-§ 

D , 


Parasite  drag,  absolute  coefficient  CDf= 


qS 


Croes-wind  force,  absolute  coefficient  C0= 


C 


& 


iw 

it 

Q 

n 

B 


a 

t 

a 0 

at 

at 

y 


Angle  of  setting  of  wings  (relative  to  thrust  line) 
Angle  of  stabilizer  setting  (relative  to  thrust 
line) 

Resultant  moment 
Resultant  angular  velocity 

Reynolds  number,  where  l  is  a  linear  dimen¬ 
sion  (e.g./foran  airfoil  of  1  .Oit  chord,  100  mph, 
standard  pressure  at  15°  C,  the  corresponding 
Reynolds  number  is  935,400;  or  for  an  airfoil 
of  1.0  m  chord,  100  mps,  the  corresponding 
Reynolds  number  is  6,865,000) 

Angle  of  attack 
Angle  of  downwash 
Angle  of  attack,  infinite  aspect  ratio 
Angle  of  attack,  induced 
Angle  of  attack,  absolute  (measured  from  zero- 
lift  position) 

Flight-path  angle 
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Page  4,  column  2,  line  3*  The  term  at  the  end  of  the  line  should 
have  a  bar  over  ^  j  thus  ?  [f>(  $1 

Page  7, 

Column  1,  equation  (28) :  The  left-hand  side  of  the  equation  should 
read  e*a 

Column  1,  line  following  oquation  (28):  Change  cr  to  9. 

Column  2,  line  16:  The  bracket  and  the  fractional  exponent  at  the 
end  of  the  equation  for  E  should  be  transposed  to  read 

...  sin  20]  1/2 

Page  13,  column  1,  last  line:  The  equation  should  read 
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ON  THE  PLANE  POTENTIAL  FLOW  PAST  A  LATTICE  OF  ARBITRARY  AIRFOILS 

By  I.  E.  Gahiiick 


SUMMARY 

The  two-tlimensional,  incompressible  potential  flow  past  a 
lattice,  oj  airfoils  of  arbitrary  shape  is  investigated  theoretically. 
The  problem  is  treated  by  usual  methods  of  conformal  mapping 
in  several  stages ,  one  stage  corresponding  to  the  mapping  of 
the  framework  of  the  arbitrary  line  lattice  and  another  signifi¬ 
cant  stage,  corresponding  to  the  Theodorsen  method  for  the. 
mapping  of  the  arbitran/  single  wing  profile  into  a  circle.  A 
particular  feature  in  the  theoretical  treatment  is  the  special 
handling  of  the  regions  at  an  infinite  distance  in  front  of  and 
behind  the  lattice.  Expressions  are  given  for  evaluation  of  the 
velocity  and  pressure  distribution  at  the  airfoil  boundary.  An 
illustrative,  numerical  example  is  included. 

INTRODUCTION 

This  paper  treats  the  problem  of  determining  the  (low  pat¬ 
tern,  or  the  velocity  and  pressure  fields,  associated  with  the 
uniform  flow  past  an  infinite  row  of  symmetrically  placed 
airfoils  of  the  same  shape.  This  airfoil-lattice  problem 
occurs  in  the  design  of  turbine  blades,  wind-tunnel  vanes  or 
grids,  and  elsewhere.  There  is  a  purely  mathematical  interest 
in  the  problem  that  concerns  the  field  of  conformal  mapping 
of  infinitely  connected  regions.  Analogous  two-dimensional 
“lattice”  problems  occur  in  the  steady  flow  of  heat  and 
electricity. 

Considerable  ingenuity  has  been  devoted  to  the  airfoil- 
lattice  problem,  especially  in  the  turbomachine  studies  in 
the  German  literature  and  more  recently  in  the  British 
studies;  nevertheless,  a  survey  of  the  available  literature 
indicates  that  nearly  all  the  treatments  employed  and  the 
results  obtained  are  of  a  special  or  indirect  nature  which 
involves,  for  example,  lattices  of  thin  lines  or  approximate 
graphical  procedures.  Recently,  however,  A.  R.  Howell  in 
a  British  paper  of  limited  circulation  has  written  briefly  on 
the  theory  of  arbitrary  airfoils  in  cascade.  Howell  applies  a 
special  transformation  to  an  airfoil  lattice  to  convert  the 
lattice  region  to  a  somewhat  random,  simply  connected 
region  and,  with  the  aid  of  several  stages  of  conformal 
mapping,  obtains  a  region  about  a  circle. 

The  problem  of  determining  the  incompressible  potential 
flow  past  an  arbitrary  single  wing  section  was  studied  by 
Theodorsen  (reference  1),  who  gave  a  practical  procedure 
for  its  solution.  The  case  of  two  wing  sections,  or  the 
arbitrary  biplane,  was  later  treated  in  reference  2.  The 
determination  of  the  flow  past  an  infinite  lattice  of  airfoils 
of  the  same  shape  is  a  problem  intermediate  in  difficulty  in 
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comparison  with  the  aforementioned  ones.  The  treatment 
for  resolving  this  problem  given  in  the  present  report  is 
similar  to  that  for  the  arbitrary  single  wing  section  but  the 
calculations  are  more  involved. 

The  problem  will  herein  be  studied  by  the  usual  method 
of  conformal  mapping.  It  is  convenient  to  accomplish  the 
result  in  three  or  four  stages:  The  airfoil  lattice  is  first  re¬ 
placed  bv  its  skeleton,  or  framework  of  line  segments.  The 
initial  mapping  function  employed  transforms  the  lattice 
skeleton  into  a  circle.  In  the  plane  of  this  circle  there  are 
two  singular  points,  known  as  branch  points.  These  points 
have  dual  significance:  They  correspond  to  infinite  regions 
in  front  of  and  behind  the  lattice  of  lines,  and  they  enter  in 
the  problem  of  reducing  the  lattice  region  (multiply  con¬ 
nected  region)  to  the  region  of  a  single  body  (simply  con¬ 
nected  region).  If  now  an  arbitrary  airfoil  shape  is  gen¬ 
erated  or  given  around  the  framework  of  lines,  then  in  the 
plane  of  the  circle  a  circular-like  contour  is  generated  around 
the  original  circle.  This  contour  may  be  transformed  into 
an  exact  circle  by  the  well-known  procedure  given  in  ref¬ 
erence  1  or  3.  The  original  two  significant  points  arc  then 
traced  by  a  transformation  due  to  H.  A.  Schwarz.  A  final 
elementary  transformation  will  bring  the  circle  into  a  stand¬ 
ard  circle  for  which  the  two  characteristic  branch  points  are 
symmetrically  placed.  The  region  of  this  circle  is  considered 
the  standard  region  for  determining  the  flow  pattern. 

For  illustrative  purposes  on  outline  of  a  procedure  for 
calculating  pressure  distributions  is  included.  The  method 
may  be  followed  without  reference  to  the  theory  by  readers 
interested  mainly  in  numerical  results.  For  convenience,  a 
list  of  symbols  is  given  in  appendix  A. 

ANALYSIS 

Initial  transformation  for  lattice  of  straight  lines. — Con¬ 
sider  the  transformation  (reference  4) 


where  g,  6,  and  a  are  real  numbers  and  bf>a.  Introduce 
coordinates  fl/  and  6  by  means  of  the  relation 

s' (2) 

and  let 


l 
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Equation  (1)  inav  than  bp  expressed  as 

-  __0_\  ,,  I" 7n  +  COsll  (<p  A/0]~I 

S1  2jt  0r-  |_cosh  To — eosli  ( ^-)-t0)  J  ' 

If  \f/= 0,  according  to  equation  (2),  lies  on  a  circle  of 
radius  a  (fig.  1(a)).  According  to  equation  (4),  u=xi-sriyi 
is  the  logarithm  of  a  real  positive  function  and,  consequently, 
represents  a  real  function  (its  principal  value)  and  the  infinite 

sequence  of  values  differing  from  this  function  bv  2kiri, 

where  k  is  any  integer.  The  transformation  illustrated  in 
figure  1(h)  is  that  of  an  infinite  lattice  of  unstaggered  lines 
of  gap  (j  in  the  fi-plane  into  the  circle  of  radius  a  in  the 
“'-plane,  'flic  points  z'—b  and  correspond  to 

infinity  in  front  of  and  behind  the  lattice,  respectively.  The 
.  .  G?  ( l * 

inverse  points  z’  —  ^  and  z'  —  ~  ■-  are  inside  the  circle  of 
radius  a. 

In  order  to  introduce  stagger,  it  is  convenient  to  consider 
the  transformation 


6-Hf 
b  —  z' 


where  h  is  real.  This  transformation  can  be  written  with  the 
use  of  equations  (2)  and  (3)  as 


■  A  t  T siQ?.1.  To + si  nh  (^  +  i0)~| 
1  2ir  og  Lsinh  To— sinh  (^-f-i0)J 


(5) 


If  ^=0,  the  expression  within  the  brackets  is  a  complex 
number  of  unit  magnitude;  hence,  the  logarithm  is  a  pure 
imaginary  number  plus  an  infinite  sequence  of  numbers 
differing  by  2iri.  Then  f2=x2-f-n/:  represents  a  sequence  of 
real  numbers  differing  by  h  and  tile  iattiee  is  one  of  hori¬ 
zontal  lines  displaced  from  each  other  by  h  (fig.  1(c)). 

The  transformation  for  the  general  staggered-line  lattice 
is  a  combination  of  equations  (4)  and  (5) 


(6a) 


or 


where 


( 


e  ^  log  log 


(Ob) 


gap  u=(l  eos  0 
stagger  h—d  sin  ;3 


stagger  ratio  — =tnn 


0 


the  parameter  d  may  he  called  the  slant  gap  (fig.  1(d)), 
and  0  the  stagger  angle. 

The  geometry  of  the  lattice  inav  he  expressed  in  terms  of 
the  parameters  y0  and  0  hv  noting  that  the  chord  length 
may  be  obtained  from  the  (singular  or  critical)  values  of  0 
which  correspond  to  the  end  points  of  the  chord  and  are 


solutions  of  the  equation  (yf,=0.  This  equation  gives  the 
result 

tan  0=tanh  y0  tan  0  (7a) 


or,  for  later  reference, 


eos  0= 


cosh  70  cos  0 


sin  0  = 


Q 

sinh  7o  sin  0 


Q 


(7b) 


where 


$=  (cosh2  70— sin2  0Y 


Relations  (7)  may  he  employed  in  two  ways:  (1)  When  the 
parameters  70  and  0  arc  given,  the  relation  determines  the 
two  critical  values  of  0,  0,  and  0,,  where  the  subscripts  l  and 
t  refer  to  leading  edge  and  trailing  edge,  respectively,  and 
0,=0i-fir.  (2)  When  0,  or  tan  0,  and  the  stagger  angle  0 

are  given,  the  relation  determines  the  parameter  y0. 


The  chord  c  may  he  obtained  by  putting  0=0 j  and  0=0, 
in  equation  (Oa)  and  taking  the  difference  in  abscissas  x,  and 
x,.  From  equations  (4)  to  (7), 


c=x,—x, 

2d  (  .  Q+c os  0  .  sin  0\  .  . 

=t  (cos  0  l0K  sinh  Tr+Sin  0  tan  ~o~)  (8) 

By  means  of  equation  (8),  the  parameter  y0  can  lie  pre¬ 
sented  directly  in  terms  of  given  values  of  the  gap-chord 
ratio  for  any  stagger  ratio.  A  representative  chart  relating 
;  gap-chord  ratio,  stagger  angle,  and  y0  is  shown  in  figure  2; 
i  some  values  are  given  in  table  I. 
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TABLE  I. — Ci AP-CIl OB D  RATIO,  PARAMETER  70,  AND  CORRESPONDING  VALUES  OF  dt  FOR  VARIOUS  STAGGER  ANGLES 


hla 

sinh  y 

a 

1*0511  y 

0 

0=0°:  tf,= o° 

j3=:io° 

die 

<i‘c 

<n  (dt’p) 

1 . 00.1 

0. 0049HH 

0.  <X)49sS 

1.  000012 

0.  20207 

0.  2948! 

0.  13 

0. 44157 

1.01 

.  009950 

.  009950 

1 .  (XNI050 

.  29019 

. 33208 

.33 

.  18904 

1.02 

.  01 9803 

.  019804 

1.  oooiuo 

.  34030 

. 37997 

.  00 

.  .14774 

1. 0.1 

. 048790 

.  048809 

1.IM11191 

.  12299 

. 10845 

1. 01 

.  0.1010 

1.  10 

.  09.1310 

.09.1455 

1.00454.1 

.  .'.1.194 

.  5004 1 

3.  14 

.  7.1.147 

1.  1.1 

.  1.19702 

.  110217 

1.009783 

.  .18995 

.  04331 

1.  .18 

. 83357 

1.20 

.  1821122 

.  183333 

1.0)0007 

,0.9.107 

.71028 

5.  94 

. 89900 

1.2.1 

.22:1144 

.  22.KXMI 

1.  02.HKXI 

.  71490 

.  77129 

7.  22 

.  9.1724 

1.30 

.  202:104 

.  20.1385 

1.  0:1401.1 

.77118 

.  82829 

8.  42 

1. 01004 

i.:i.i 

.  300105 

.  304030 

1.  (145370 

.  >2489 

>8241 

9.  .15 

1.00073 

1.40 

,  330472 

. 342857 

1.057113 

.  87008 

.  934:i3 

11).  01 

I. 10841 

1.45 

.  .171504 

. 380172 

1.009828 

.  92094 

.  984.10 

1 1 . 59 

1. 15433 

1.50 

.  40.1405 

. 110007 

l.08:i:i:i:i 

.  97599 

1.03340 

12.52 

1. 19890 

1.00 

.  171X104 

.  IS7.KX) 

l.  112.KX1 

1.07124 

I. 12792 

14.20 

1.  28505 

1.70 

.  &KM2H 

.  .1.1.1882 

I  1-14  UK 

1. 10302 

1.  21927 

1.1.  07 

1. 30849 

1.  HO 

.  .1*7787 

1.  177778 

1.2.1387 

1 . 30830 

10.  90 

1.  4.1020 

1.  DO 

.•►41854 

.  080842 

1. 213158 

1.342.10 

1,  39501 

IS.  10 

1.  53007 

2.  0 

.  09:1147 

.  7.KXXXI 

1 . 2.KXXX1 

1. 42980 

1.48157 

19.  11 

1.01030 

.  910291 

I.  O.KNNXI 

1. 1.KXXXI 

1.85391 

1.  s98m 

22. 09 

2.  (X)292 

.Lo 

i.0980i:t 

1.333333 

1 . 0000117 

2.  20019 

2. :  10.153 

24.  79 

2.  39240 

4.0 

1.  .180294 

1 . 87.HXX) 

2.  12.9000 

3.  07484 

3.  10.188 

27.  (XI 

3. 171  IS 

.1.0 
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Inversion  of  equations  (4)  to  (6). — Tho  initial  transforma¬ 
tions  may  ho  thought  of  as  mapping  a  framework  of  chords 
of  an  arbitrary  lattice  into  a  circle.  If  a  shape  is  generated 
around  the  chords  in  tho  z'-plano,  a  contour  is  generated 
around  the  circle  of  radius  a.  This  contour,  which  must 
exclude  the  points  z'  —  —  b  and  2' =  6  and  must  enclose  the 
at  ft2 

points  z'=  — -j-  and  2'  =  ^>  >na.v  ho  considered  to  he  com¬ 
pletely  defined  by  the  function  4/(8).  If  a  lattice  of  airfoils 
is  proassigned,  the  function  4/ (9)  must  he  found  from  tho  given 
coordinates  of  the  airfoil  shape.  In  order  not  to  interrupt 
tho  sequence  of  main  ideas,  the  details  of  this  problem  are 
relegated  to  appendix  B,  with  certain  remarks  on  the  practical 
achievement  of  a  nearly  circular  contour. 

Transformation  of  contour  in  s'-plane  to  circle  in  2- 
plane. — It  is  assumed  now  that  the  circular-liko  contour 
in  tho  s'-plane  which  corresponds  to  the  airfoil  contour  of 
the  lattice  is  either  given  or  determined ;  that  is,  the  function 
4/(8)  is  known  in  the  boundary  expression  z'=ae*+,e.  By 
the  procedure  of  reference  1  or  3,  the  transformation 

2'  =  zefiz)  (On) 

where 

./(2)=S!»= 102:7-  (9b) 

and  c„  are  complex  coefficients  determined  by  tho  boundary, 
is  then  employed  to  transform  tho  V -contour  into  a  circle 
in  the  2-plane.  The  transformation  (Oa)  keeps 
the  regions  alike  at  infinity  in  the  z'-  and  r-planes;  that  is, 

(lzr 

2=2'  and  ^-=1  at  infinity.  The  correspondence  of  the 
boundaries  is  determined  by  (lie  functional  equation 
4>—8=t(<t>) 

=  d<j>'  (10) 

for  which  a  convenient  numerical  solution  has  been  outlined 
in  reference  3.  The  radius  of  tho  circle  ll=ae+a  is  determined 
by  the  relation 

1  ... 

I  (if) 

-hr  J  i) 


For  consistency,  the  functional  symbol  ^(1 j>)  is  here  used  to 
denote  the  quantity  4'  expressed  as  a  function  of  4>  —  that 
is,  1 p[8(4>)].  In  reference  3  the  notations  4/{4>)  and  i/-[0(<6)] 
are  used. 

It  is  necessary  also  to  trace  the  correspondence  of  the 
points  z'=b  and  s’  =  ~b.  Let  2=0 ,  correspond  to  z'=b 
and  let  2=  —  02  correspond  to  z'——  b.  The  values  0,  and 
0o  may  he  determined  by  a  relation  (due  to  Schwarz)  that 
expresses  the  value  of  a  complex  function  in  terms  of  an 
integral  of  the  real  part  of  the  function  along  a  circle.  A 
simple  derivation  of  the  desired  relation  is  shown  in  appendix 
C.  The  expression  is 

log  j=.f(z) 


The  values  of  0,  and  02  may  he  determined  from  equation 
i  (12)  hv  an  iteration  process  that  converges  extremely 
I  rapidly.  The  process  may  be  described  as  follows:  In 
equation  (12),  let  tho  zeroth  approximation  to  0i  be  z=z„=b 
and  let  the  corresponding  value  of  z’  be  written  z'  =  20' =befm , 
where  /(6)  is  the  evaluation  of  equation  (12)  for  2=6.  It  is 
actually  desired,  however,  to  have  z’=b  hut,  because 

z'  =  20'  =6+  2u'  —  6 

tho  initial  value  of  2'  differs  from  the  desired  value  l>v 
2,/  — 6.  Furthermore,  2=r0  differs  from  2=0 x  by  approxi¬ 
mately  the  same  amount;  hence,  reducing  20  by  the  quantity 
Zj—b  gives 

Zi=z0+b— 2o' 

=6[2— e/(M] 

which  may  he  considered  a  first  approximation  to  0,.  If  it 
is  desired  to  check  this  result  or  to  obtain  a  second  approxi¬ 
mation,  the  process  may  he  repeated;  thus,  from  equation 
(12),  find  j(zx)  and 

2, '  =  2,^'' 
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which  is  a  second  approximation  to  ft  and,  in  general,  the 
a th  approximation  is 

n  —  ‘•n—  1  d~  5  — n— 1 

It  is  clear  that,  should  z„  correspond  to  r„_/  must  cor¬ 
respond  to  6  and  the  process  automatically  stops.  The  nu¬ 
merical  process  is  given  in  appendix  O;  relatively  elemen¬ 
tary  steps  are  involved.  In  order  to  determine  — ft,,  the 
process  is  applied  with  6  replaced  by  —h. 

Transformation  to  standard  circle  in  it- plane. — In  order 
to  obtain  the  llow  pattern,  it  is  desirable  to  introduce  another 
function  which  transforms  the  circle  in  the  r-plane  into 
another  circle  in  the  « ’-plane  in  such  a  way  that  the  char¬ 
acteristic  i>oints  z=\ 3,  and  2=  —  ft  map  into  w—b  and  to=  —b, 
respectively.  The  region  of  the  circle  in  the  tr-plane  mav 
be  considered  the  standard  region.  The  desired  trans¬ 
formation  may  be  written  ns  (see  appendix  D) 


b~w  k(01~2\ 

b+w  \ft+2/ 

where 

rsb^2  ftft  +1? 

ftft~/r: 


(13) 

(14) 


and  R=ae*t 1  is  the  radius  of  the  original  circle  in  the  2-plane, 
ft  is  the  complex  conjugate  to  ft,  and  S  is  the  radius  of  the 
new  circle  in  the  w-plane.  The  radius  S  is  determined  bv 


lie  readily  verified  that  the  circle  of  radius  — that  is, 
w=Sei'’ — is  part  of  a  streamline  and  it  may  further  be 
observed  from  figure  3  that  the  circulation  around  any 

S2 

contour  which  encloses  the  points  w=  ±  and  for  which 

the  points  w=+b  are  exterior  points  is  f  (positive  if  coun¬ 
terclockwise).  The  parameter  a  will  be  interpreted  later 
as  an  angle  of  attack. 

The  value  of  the  circulation  T  may  ho  determined  by 
means  of  the  Kutta-Joukowski  condition  for  smooth  llow 
at  the  trailing  edge  of  the  lattice.  Let  aa  be  the  value  of  a 
on  the  boundary  circle  Se*”  that  corresponds  to  the  trailing 
edge  of  the  lattice.  The  Kutta-Joukowski  condition  then 
requires  that  the  flow  separate  at  o~an,  or  that  a  stagnation 
point  exist  there. 

(IQ 

With  (j'w= (•  and  w—Se11”1,  the  following  relation  for  f  is 
found: 

r=—  •— yf[^  sin  (<r°-l-a)-(-y  sin  (<r0— a)J  (IS) 

h'~  )? 


If  Sib  is  replaced  by  e~y'  (equation  (1*0),  e<|uation  (IS)  may 
be  expressed  as 


S=6«-r' 


where  y,  is  obtained  from 


cosh  7! = 


ll 


fl2+ftft 
ft  +  ft 


do) 

(Hi) 


Complex  velocity  potential  in  w-plane. — Consider  the  flow 
function  V.(w)  =<l>+t'l',  which  is  defined  as 


<-;(«•) 


log 


b+w 

b—\o 


+e*‘“ 


log 


b‘ —  MT 


(17) 


'flic  flow  pattern  may  be  regarded  as  due  to  a  combination 
of  singularities,  sinks,  sources,  and  vortices,  placed  at  the 
S2 

points  w=  ±b  and  w=  ±  y  as  indicated  in  figure  3.  It  may 


Expressions  for  velocity  in  lattice  field. — Tn  order  to  obtain 
the  flow  pattern  in  the  lattice  field  (f-plane),  the  component 
factors  of  the  following  expression  are  required: 

(Kl_<lndw<lz  dz!_  ,0f.. 

d{~dw  dz  dz'  dt;  ' 

These  terms  may  be  obtained  from  equations  (17),  (13), 
(9),  and  (0). 

It  is  of  particular  interest  to  evaluate  equation  (20) 
explicitly  for  the  regions  at  infinity  in  front  of  and  behind 
the  lattice  and  also  on  the  lattice  boundary  itself.  It  is 
recalled  that  corresponds  to  z'—b,  z  =  /So  vs=b  and  that 
f=  —  co  corresponds  to  z'  —  —b,  z——  ft,  w=—b.  By  com¬ 
bining  terms  according  to  equation  (20),  the  (reflected) 
inlet-velocity  vector  is  obtained  as 


iT 


=  _  IV<a+tfi  —  — ,  e. 

-Id 


(21) 


and  the  corresponding  expression  for  the  outlet-velocity 
vector  s 


W1 

di  J- 


= V„-iV 


yTgi(a+8> 


2d 


(22) 
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By  mldiiioti  of  equations  (21)  and  (22),  il  becomes  clear 
tlint  l ho  velocity  vector  of  magnitude  l"  and  angle  of  attack 
a  +  /3  with  respect  to  the  /-axis  is  one-half  tile  vector  sum 
of  tiie  inlet  and  outlet  velocities  (fig.  4). 


Fioure  4.— Inlet,  outlet,  nml  moan  velocity  vectors  and  angles  of  attack. 


If  the  angle  of  attack  of  t tie  mean  velocity  vector  with 
respect  to  the  g-axis  (chord  direction)  is  denoted  by 
az=a- \-0,  the  velocity  components  in  equations  (21) 
and  (22)  arc 


VX{  =  —  V  cos  az  +^  sin  6 

r 

Vvx  =  —  V  sin  az  cos  0 


by  angle  8  (lig.  4).  These  components  are.  for  the  inlet 
velocity, 

1  Vi  —  —  1  COS  Cl 


l 


£,=  V  sin  a  — 


2d 


and,  for  the  outlet  velocity, 


I  ,v„=  —  I  COS  0—1  y  [ 


rl2=  I'sin  a~—j 


The  squares  of  the  magnitudes  of  the  inlet  and  outlet 
velocities  are 

[' +22!a  «-(_.!•,/)  ] 

where  17217/  may  he  obtained  from  equation  (19).  Observe 
that  the  inlet  and  outlet  speeds  are  equal,  r,=  I'2>  when 
a=0°  for  any  value  of  T.  The  inlet  and  outlet  angles  of 
attack  with  respect  to  the  normal  to  the  lattice  line  arc 


m=tan-‘  sin  a+Wd 
cos  a 

V 

■  sin.  ot  t  » » 
a2=tan  1 _ 2  V a 

cos  a 

and  the  angle  through  which  the  stream  is  turned  is 


c*!  —  a2  =  tan  1 


o 


2lWCOSa 


(23) 


The  component  factors  in  equation  (20)  are  now  to  be 
evaluated  at  the  lattice  boundary  and,  as  the  boundary 
itself  is  part  of  a  streamline,  only  the  magnitudes  of  the 
factors  are  of  interest. 

From  equations  (17)  and  (19)  and  with  w=Se‘1’, 


dV 

ilw 


2 Vd _ 1 

T&  0081127,-003  2(7 


[sinh  7,  sin  a(eos  a— cos  <t„) 


and 


i'V.=  -r 


COS  ar~.,^  sill 


ly,,  =  V  sin  az  —  :^i  cos  8 

The  conventional  angle  of  attack  a  is  measured  with 
respect  to  the  normal  to  the  slant  line  of  the  lattice.  The 
components  normal  to  and  along  the  slant  line  of  the  lattice, 
sometimes  referred  to  as  “axial”  and  “whirl’’  components, 
respectively,  are  found  by  rotating  all  v<  ctr.rs  in  the  /’(/-plane 


+  cosh  7,  cos  a(sin  o—  sin  <r„)]  (24) 


where  the  parameter  7,  is  defined  in  equation  (15). 

In  order  to  obtain  dwUlz,  it  is  convenient  first  to  express 
equation  (13)  explicitly  in  w  as 


6(1+  K)  z  -6  (Kg, -ft) 

(1  —  A*)2-t-A'ft  +  ft 


(25a' 


A  standard  form  for  the  transformation  of  a  circular  region 
;r  S  R  into  |roj  7  b’  is 


w=HSe» 


:— 8_ 
H-—h 


(25b) 
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Comparison  of  equations  (25a)  and  (25b)  makes  it  clear  that 
the  complex  parameter  5  and  the  real  parameter  X  may  he 
obtained  from  the  following  relations: 


,  A'ft-ft 

&-~T+K~ 


(20a) 


or,  as  a  cheek  relation, 

.  R2(K-l) 

6  A'ft  +  ft 

and 

S  n  ( 1  + 10  b 

Ile  ~A'ft+ft 

or,  by  equaling  angles  on  both  sides, 

X = a  rg  ( 1  +  A”)  —  a  rg  ( Aft + ft) 


(20b) 


(27) 


From  equation  (25b),  the  explicit  correspondence  of  a  point 
on  the  eirele  w=Sef’  to  a  point  on  the  circle  z=Rei*  can  be 
obtained  as  follows: 


pi  ($+Xi_ 


Re 

I 

1  Re 


(28) 


Let  the  complex  number  a  be  expressed  as  \S'eir  and  let 


where 


e~i*=meil‘ 


jw  '^,2 

=  1—2  ^  cos  (</>  — t)+^2 


(29) 


and 


sin  (<t>—  r) 

nW  =  tan-1 - ttt - 

1— • ^  cos  (</>— r) 


Observe  that  the  denominator  in  equation  (2S)  is  the  con¬ 
jugate  of  equation  (29)  and  is  therefore  equal  to  me'1*. 
There  results  for  the  correspondence  of  <r  and  </> 

v=<f>-(~  X+2g  (30) 


In  particular,  if  the  (trading-edge)  value  of  <t>  that  corre¬ 
sponds  to  6,  as  determined  by  equations  (7)  is  written  as 
4>0—e,+ti,  where  e,  is  the  value  of  «(</>)  at  6=6,  from  equation 
(10),  then 

<70  =  <f>0  +  ^  +  2/xo 


By  differentiation  of  equation  (256), 

dw  _  RS(R2—sd)e'>' 
dz  (A2— js)2 


CU) 


On  the  boundary,  put  z=Reu;  then,  the  magnitude  of 
equation  (.31)  is 


dw._  S  (  _j5/\  J_ 

■dzrRY  R2)  m~ 


(32) 


rir.aTi  — 4S - 2 


The  expression  for 


dz"' 

dz, 


on  the  boundary  is  obtained  from 


equation  tit)  in  terms  of  the  functions  ei</>)  and  'P(</>)  of 
equation  (10)  as  follows  (see  reference  3): 


dz' 

dz 


(33a) 


and,  because  f(z)  on  the  boundary  is 

/(c)  ='*'(<#’)  —  <l/o+i(6—(t>) 


where 

then 


'dz'  | 

uizY 


z  [(>4;HS)T 


(33b) 


The  last  factor  of  equation  (20)  is  expressed  from  equation 
((>)  on  the  boundary  s'=ae*+M  as 


where 


</f  =2d  Al_ 
dz'  r  U  z' 


(34) 


A=|^cos2/3  cosh2y0(cosh2^— cos20) 

1  Hi 

-r  sin2£  sinh27o(cosh2i/'— sin20)  — —  sin  2/3  sinh  2y0  sin  26 

D=c osh  2y0—  cosh  2(\f/+i6)\ 

=  [(cosh  27o— cosh  2 cos  20)2+  (sinh  2^  sin  2 0)2]t* 

Finally,  combining  in  equation  (21)  the  factors  given  in 
equations  (24),  (32),  (33b),  and  (34)  yields 


=  ABCD  V  (35) 

where 


-1  =  —  ,—7: - x-  [sinh  7.  sin  a  (cos  a — cos  vo) 

cosh  27,  —  cos  2  cr  1  u/ 

—  cosh  71  cos  a(sin  a—  sin  <r0)] 

c-[(‘-/;)’+('S)T‘ 

0=1  (cosh  270— eosli  2\p  cos  20)*+ (sinh  2i p  sin  20)2]-4 
A= J^cos2/3  cosli27o  (cosh2^ — eos20) 

1  .  *l>s 

-rsin-’d  sinh27i,(coslr'i/'— siri20)  —  —  sin  2/3  sinh  270  sin  20 

An  application  of  equation  (35)  for  the  purpose  of  illus¬ 
trating  the  various  steps  involved  in  a  calculation  of  the 
surface  velocity  and  pressure  of  the  airfoil  lattice  is  given  in 
appendix  E  and  illustrated  in  figures  5  and  0.  For  the  sake 
of  comparison,  the  single-airfoil  case  is  given  in  figure  7. 


s 
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Floi'RR  5.—  Pressure  distribution  for  N  A C  A  -4412  airfoil  in  a  lattice  arrangement.  Stunner  uncle  0-0°:  cap-chord  ratio,  O.UiW;  slope  of  lift  curve 
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Some  special  results  from  equation  (35)  for  a  lattice  of 
lines. — Til  the  case  of  a  lattico  of  straight  lines,  the  z'-, 
and  w-planes  merge;  lienee  0=<t>—<r  and  R=S=a. 

From  equations  (10)  and  (7)  and  with-  a+0=ax,  which  is 
the  angle  of  attack  with  respect  to  the  chord, 


2Vd  (cosh2  70— sin*  0)^  2 

The  lift  per  unit  span  on  a  single  member  of  the  lattice  is 
given  by 

L=pVV 

where  p  is  the  air  density.  The  lift  vector  is  perpendicular 
to  the  mean  velocity  vector  (fig.  4).  This  result  is  general 
and  not  limited  to  a  straight-line  lattice.  The  lift  coeffi¬ 
cient  is 


„  P^r 


<>■) 


_ 2r_  _i _ r 

cV  cjd  21  d 


_ 1 _ .  cos  0  eos  9  .  sin  0  sin  9 

(cosh2  To  —  sin2  0)H cosh  70  sinh  y0 

^  eos  0  sin  0  sin  0  eos  6 
1  sinh  70  cosh  70 

In  the  special  cases  in  which  0=0°  and  /3=90°,  the  relations 
(30)  to  (38)  arc  simpler. 

For  stagger  angle  0=0°  and  with  d=(j, 


From  equation  (8), 


T  _  sin  aT 
2  Vg  cosh  7o 


cosh  7n=coth  * 


L=2pV1g  tank  ~  sin  ax 
tanli  ™ 

=  ir pcV2 - —  sin  a, 

ire 


where  T/2  Vd  is  given  in  equation  (30)  and  c /d  can  be  found 
by  equation  (8). 

The  local  velocity  on  the  surface  (equation  (35))  becomes 


»=V^cos  <*r-|-^sin 


The  lift  coefficient,  according  to  equation  (37),  is 

tanhrp 

CL=2ir - —  sin  ax 

JT6‘ 


1U 


REPORT  No.  7S8 — N  ATIONAL  ADVISORY  CUlt.M  ITTKK  FOR  AERONAUTICS 


For  ii  —  0°,  therefore,  the  slope  of  the  lift  curve  is  always  less 
than  2v.  Note  that,  for  la  [fro  yap,  <\p— 0  and  the  lift 
eoellieient  is 

Cl='2tt  sin  a. 

When  the  gap  <j  is  small  compared  with  the  chord  c, 

ii  . 

(  t  sin  ctT 

The  local  velocity  at  the  surface,  by  equation  (38),  is 

v—  Pr^cos  aT+tnnh  7,1  cot  ~  sin  a,j 

This  result  may  he  compared  with  that  for  the  single-line 
airfoil  (7,)=“) 

v=  r(cos  a, -foot  7  sin  a, ^ 

For  stagger  angle  /3=f)0°  and  with  d=h, 

F  sin  a. 


From  equation  (8), 


2  Vh  sinh  7j 


sinh  y0=coty| 


and 


L  =  2pV2h  tan  sin  a, 

,  TC 

l/2  an  2h  . 

-  Trpc  V  2  - sin  az 

7 rc 

•7* 


The  lift  coefficient,  according  to  equation  (37)  is 


4  7TC 

r  „  tan2A  . 

(,V=2x  - stn  a, 

vc 

2h 


i  Foi  [i - -00°,  therefore,  the  slope  of  the  lift  curve  is  always 
j  greater  than  27r.  The  local  velocity  at  the  surface  is 

r  =  I’jj'os  «,-j-coth  7 ,,  cot  ^)sin  arJ 


It  may  he  noted  in  passing  that,  for  c  =  y/i , 
CL=S  sin  a. 


as  compared  with 


(\=2ir  sin  a. 


for  the  single  airfoil. 

For  the  limiting  case  in  which  b  and  d  approach  «> ,  the 
transformation  (0)  becomes 


-=  d  al  \ 
•  2  irb\ 


and,  with  limit  , »1  and  a  new  variable  2"  =  s'e~is 
2vb 

f=z*+5 

* 

which  is  the  familiar  Joukowski  transformation.  If  the 
variables  <!/  and  8  (equation  (2))  arc  introduced,  the  cor¬ 
responding  result  is  expressed  as 

f=2 a  cosh  lt+i(0— 0)] 

where  the  limit,  as  y0~*°° ,  of  y - —  has  been  put 

’  ,u  ’  2irO  cosli  70  1 

equal  to  1. 
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APPENDIX  A 


MAIN  SYMBOLS 


if  complex  plane  of  airfoil  lattice  (x  +  iy) 

fi,  {".,  complex  planes  of  airfoil  lattice  for  stagger  angles 
(3  =  0°  and  /3  =  9<)°,  respectively  (jvt-d/r,  ft+i'h) 
■  complex  plane  of  circular-like  contour  (ae*+ii) 

:  circle  of  radius  R—ue*'1  in  z-plane 

w  circle  of  radius  S—be~yt  in  w-plane  (be~y'+,‘’) 

f=ao,  z'  =  b,  c=(3i,  w = h  corresponding  points 

f=  — oo,  z'  =  —  b,  z—  —  i 3-2,  w—  —  h  corresponding  points 

«.  b  reference  lengths 

7o  gap-chord  parameter  (b—oe1  o) 

i3  stagger  angle 

il  lattice  spacing,  or  “slant”  gap  for  any  value  of  d 


y  lattice  spacing,  or  gap  for  (3  =  0° 

It  lattice  spacing,  or  stagger  for  (3  =  90° 

V  magnitude  of  mean  of  inlet-  and  outlet-velocity 

vectors  (jig.  4) 

a,  angle  of  attack  with  respect  to  /-axis  of  mean 
velocity  vector 

a  angle,  of  attack  with  respect  to  normal  to  slant  line 

of  lattice  of  mean  velocity  vector 
on,  a2  inlet  and  outlet  angles  of  attack  with  respect  to 
normal  to  slant  line  of  lattice  respectively 
V,,  1  j  magnitudes  of  inlet  and  outlet  velocities,  re¬ 
spectively 


APPENDIX  B 


INVERSION  OF  EQUATIONS  (4)  TO  (6)  AND  CHOICE  OF  COORDINATES 


It  is  desired  to  find  from  a  given  airfoil  lattice  in  the 
{•-plane  the  contour  defined  by  \p( 6 )  in  the  z'-plane.  This 
problem  corresponds  to  an  inversion  of  equations  (4)  to  (0) 
and  can  be  exactly  treated  for  the  cases  in  which  (3—0°  and 
(3=90°  (equations  (4)  and  (5),  respectively)  but  an  iteration 
or  successive-approximation  method  is  required  for  equation 
(G).  Furthermore,  although  the  parameters  g  and  h  are 
fixed  by  the  geometry  of  the  lattice,  a  choice  exists  in  the 
definition  of  the  chords  and  the  origin  of  coordinates.  This 
choice  is  discussed  following  equation  (B17). 

Stagger  angle  (3  =  0°. — From  equation  (3),  there  is  obtained 


cosh  (^  +  i0)=cosh  7o  tunh  -  {4 


(Bl) 


Putting  {■[=/,+()/!  and  denoting  the  real  and  imaginary 
parts  of  equation  (Bl)  by  and  y, ,  respectively,  leads  to 


sh  4’ 


cos  #=£,= 


cosh  7o  sinh  ~T  x, 


cosh  7t>  sin  -  -  iq 
sinh  4*  sin  0=%  = - r, - — 


cosh  —  z,  4-eos  -  -  >/, 
!)  U  ’ 


( 132) 


'I'lie  expressions  containing  z,  and  1/1  in  equation  (B2)  are 
considered  given  since  these  quantities  are  known  from  the 
coordinates  of  the  airfoil  lattice.  If  41  and  d  are  eliminated 
successively, 


and 


■(-&-)  =  1 
\sin  e) 


( -li — V- 

\cos  8/ 

\cosh  \p/  ”\sinh  4'/ 


(B3) 


From  equation  (B3),  there  result  the  following  expressions, 
which  serve  to  define  the  function  \p(9)  in  terms  of  tho  airfoil 
coordinates: 


sin2  8=2>+\i>i+ru1  | 
sinh2  \!'  =  —  P+vV-HfI 


(B4) 


where 


V 


1 

•» 


(l-fi’-ij.2) 


For  small  values  of  9,  the  relation  sinh  mav  be  used. 

’  sin  8  J 

It  is  useful  for  computational  purposes  to  record  the  real 

and  imaginary  parts  of  equation  (3) 


when* 


?/i=^  (<3>i  —  <fe)  j 


(Bo) 


p2—  (cosh  70+cosh  cos  0)2-j-(sinh  <P  sin  8)2 
p,2=(cosli  70— cosh  cos  0)2+(sinh  >£  sin  0)2 
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sin  4>i=  sinh  0  sin  8 
Pi 

sin  0>  =  — -  sinii  0  sin  0 

P 2 

The  angles  are  to  he  chosen  between  —  x  anil  x,  anil  the 
quadrants  may  be  determined  by  noting  also  the  relations 

cos  <i>i=  -  (cosh  7n+eosh  0  cos  6) 

Pi 

cos  0i=  -  (cosh  y0— cosh  0  cos  6) 

pj 

Stagger  angle  /3=90°.— From  equation  (5),  there  is 
obtained 

sinh  (0-H0)=sinh  ylt  tan  j  f2  (BC) 


With  fa=J\>+i?/ 2  and  the  real  and  imaginary  parts  of  equation 
(B6)  denoted  by  £>  and  ri2,  respectively, 


sinh  0  cos  0=£.  =  - 


sinh  70  sin  ^  x2 


cosh  V  y2+cos  V  x2 


sinh  7o  sinh  y. 


cosh  0  sin  fl=ij2  = 

cosh  ^  y2+cos  go 
If  0  and  0  are  eliminated  successively, 


(Jl 

\cos  0 


*_y 

\sinh  0/  \eosli  0/ 


(B7) 


(B8) 


From  equations  (B8)  there  result  finally  the  following 
expressions,  which  serve  to  define  the  function  0(0)  in  terms 
of  the  airfoil  coordinates: 


where 


cos20=y+ vV+fc2 

sinh20=—  y+  -\V+£.'2 


(l-fc2-*2) 


(B9) 


For  values  of  0  near  ±90°,  the  relation  sinh  0=  J^2-  mav 

cos  0 

be  used. 

It  is  useful  for  computational  purposes  to  write  the  real 
and  imaginary  parts  of  equation  (5) 


X2  (03  —  0l) 

A/1,  P32\ 

=  2x  \  2  °”  p?/ 


Ih 


(BIO) 


Pa'=  (sinh  7,1-rsinh  0  cos  0)2-f-  (cosli'0  sin  0)2 
p,-=  (sinii  7„— sinh  0  cos  0)“+  (cosh  0  sin  0)* 


sin  0:,=  -  cosh  0  sin  t 

P3 


sin  0,  =  —  cosh  0  sin  0 
Pt 

The  angles  are  to  be  chosen  between  — x  and  x,  and  the 
quadrants  may  be  determined  by  noting  also  the  relations 


cos  03=  (sinh  7,,-f-sinh  0  cos  0) 

P3 


cos  0i=  (sinh  70  — sinh  0  cos  0) 

Pi 

Arbitrary  stagger  angle  0  and  choice  of  coordinates. — 

Because  of  the  transcendental  nature  of  equation  (6),  a 
direct  inversion  expression  seems  unobtainable;  however, 
the  values  (0,  0)  that  correspond  to  coordinates  (x,  y)  may 
be  obtained  without  difficulty  by  an  iterative  process.  For 
this  purpose  and  for  the  purpose  of  choosing  the  coordinate 
axes,  expansions  of  x,,  x2,  ?/,,  and  i/2  in  powers  of  0  are  useful. 
The  following  expansions  may  be  readily  verified: 


Xi  cos  0 


[M 


osh  7n+ cos  0 
osh  7„ — cos  0 


,,2  i  „  sinh2  7„— sin2  0  , 

+  02  cosh  7o  COS  0  7 - f- - rHviT-  .  .  . 

(coslr  Yo— COS2  0)2  J 


X2 


«  --  sin  0  2 


tan’ 


sin  0 
sinh  Yo 


•  „  cosh2  7,,-rcos2  0  , 

+  02  sinh  70  sin  0  - - x-tt-.H-  .  .  . 

(cosh-  7u — cos-  0)2 


?/i  «  —  cos  0 


Then 


where 


2  cosh  7n  sin  0 
cosh2  70— cos2  0  ^ 

<1  ■  „  2  sinh  7o  cos  0  , 

Vi  ~  ~rr  sm  d  — ,  . - t~o  0 

■'  2  x  cosh-  Yo— cos2  0 

y—in+ih 

0F(0) 


(Blla) 

(Bl  lb) 
(Bile) 
(Bl  Id) 

(B12) 


*’(*)  = 


cosh  7n  cos  d  sin  0— sinh  y,,  sin  d  cos  0 
cosh2  Yo— cos2  0 


If  the  x-coordinate  of  the  straight-line  lattice,  which  is  con¬ 
sidered  the  skeleton  of  the  airfoil  lattice,  is  denoted  by  Xq, 
then  x0  is  given  by  the  value  of  x=x1+x2  for  0=0,  or 


d  (  a  ,  cosh  Yn+cos  0  ,  sin  0  \ 

Xo=,t-(  cos  0  log  — j - -^+2  sin  0  tan  1  —r —  ) 

2x\  b  cosh  yo— cos  0  sinh  Yo/ 


and 


***+£  f-G(6) 


(B13) 

(Bid) 
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where 


With  this  value  of  ?/„,  and  equation  (Bio), 


0(0)- 


oosli  7„  cos  0  cos  04-sinh  jo  sin  0  sin  6 
coslr  70— cos2  6 


_  (cosh  7n  cos  0  sin  6— sin  y„  sin  0  cos  6)  2  sin  6  cos  6 
(cosh2  7o— cos2  0)2 

=  t’W) 

In  particular,  the  leading-  and  trading-edge  points  x=x,  and 
x=x,,  arc  determined  by  the  values  of  6=6,  and  6  =  6,  that 
may  be  obtained  from  equations  (7b).  Then, 


x—So, 

If  the  total  ordinate  for  both  upper  and  lower  sides  at 
r.=Xvt  is  denoted  by  ?/,, 


\h 

Xl—X o, 


*4 


(BIT) 


where 


X I  “  4'Gn 


G0=Q  (-f0?!  “  +  ~  c> S  ) 

c  \smh-  7„  cosh*  yj 


(B15) 


and  x0,  denotes  the  leading  edge  of  the  line  given  by  ^=0. 
A  similar  expression  holds  for  x,. 

From  equation  (B12),  for  constant  i p, 

-if  am 

In  the  neighborhood  of  the  leading  edge,  therefore, 

y^G0(6-6,)  (B16) 

7T 

For  x0  near  x0,,  there  is  obtained  from  equation  (Bl3), 

Xo=*ol+  {f,—6i)x 0|'  +  -f1'*  x„"  +  .  .  . 

where  the  following  relations  are  found  to  hold: 


[*■(*)]  «-*,= o 

Xo  ,"--[-0(0)],., 

7T  ‘  7T 


Hence. 


Xu~Xo,  G>,(6 —  6  ,)2 

Then,  from  equation  (B 14), 

x— x0,*x0— x0,+or  f'G(6) 


«^0„[-(0-O,)2+^l 


It  follows  from  equation  (BIG)  that.  forx=xCii, 

0-0,z=\P 

and 


This  result  leads  to  a  simple  and  convenient  way  of  choosing 
axes  of  coordinates  in  order  that  4/(6)  will  behave  smoothly 
at  the  edges,  that  is,  that  the  value  of  4/  at  the  leading  edge 
is  approximately  the  mean  of  the  values  of  4/  at  nearby 
ordinates  on  the  upper  and  lower  surfaces.  For  a  parabola 
the  latus  rectum,  or  ordinate  through  the  focus,  is  four  times 
the  distance  from  the  vertex  to  the  focus.  Equation  (B 17) 
states  that  the  end  point  of  the  skeleton  chord  should  be 
approximately  the  focus  of  a  parabola  at  the  nose. 

The  scheme  for  choice  of  axes  is  as  follows:  Locate  a  point 
Fnoar  the  leading  edge  where  the  ordinate  through  Fis  four 
times  the  distance  of  F  from  the  leading  edge.  Similarly 
locate  a  point  F'  near  the  trailing  edge.  The  origin  of  coor¬ 
dinates  then  bisects  the  line  FF' ,  which  is  on  the  x-axis  and 
represents  the  chord  of  the  skeleton  line  airfoil  4>— 0.  (To 
the  order  of  approximation  employed,  the  aforementioned 
choice  of  axes  coincides  with  that  given  for  the  single  wing 
section  in  reference  1  or  3.) 

Procedure  for  finding  (4,  6)  from  (x,  ?/)  for  arbitrary 
stagger  angle  0. — An  iterative  procedure  is  given  herein  for 
finding  4/(6 )  from  (x,  y)  for  arbitrary  0,  in  which  the  knowledge 
of  the  case  for  0=0°  is  employed  to  help  in  formulating  the 
initial  approximation.  In  brief,  values  of  6  are  obtained  for 
stagger  angle  0=0°  for  both  the  airfoil  and  its  line  skeleton. 
Values  of  6  are  then  found  for  the  skeleton,  in  the  case  of 
stagger  angle  0.  These  functions  permit  approximate  values 
of  6  to  be  found  for  the  airfoil,  for  stagger  angle  0.  Equation 
(B12)  then  enables  approximate  values  of  4  to  be  obtained. 
These  values  of  (4,  6)  are  then  readily  checked  and  improved, 
if  necessary.  The  steps  arc  as  follows: 

(1)  Choose  the  axes  as  outlined  and  express  the  airfoil 
coordinates  in  percent  chord,  where  the  chord  for  this  pur¬ 
pose,  is  the  part  of  the  x-axis  intercepted  by  the  airfoil. 
Denote,  the  coordinates  thus  obtained  by  (x„,  yp).  Find 
k=FF'  in  percent  chord.  Find  x,—xa„  the  distance  from  the 
leading  edge  to  F  in  percent  chord,  and  denote  this  value  by  e. 
Obtain  the  ratio  c/d,  where  c  means  here  FF'  and  d  is  the 
spacing  between  corresponding  points  on  adjacent  airfoils  of 
the  lattice.  Find  conversion  factor  m  by 
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(2)  Convert,  coordinate's  of  tin*  airfoil  from  {/.„  //,,)  to 
2 7 r  jji  2 ir  ^  ^  as  follows: 

k 


2 7T  M 


(e+r,-+>) 


..  y 

(3)  Find  t ho  parameter  y„  that  corresponds  to  the  deter¬ 
mined  value  of  c/d  for  the  given  value  of  0  from  graph  or  by 
ealeulation  (equation  (8)).  Also  find  for  later  use  the  value 
of  c/y  corresponding  to  this  value  of  yn  for 

(4)  Consider,  for  this  value  of  yn,  the  two  straight-line 
eases  (0=0,  0=0°)  and  (0=0,  0=0)',  associate  values  of 
0=0, ,  for  0=0°  with  values  0=0 n  for  the  stagger  angle  0  by 
referring  associated  values  of  0  to  geometrically  similar 
[joints  of  the  lines  (equation  (Bid)). 


(5)  Multiply  coordinates  in  step  (2)  by  the  ratio 


£<;/>/)  o 

where  the  chord-gap  values  are  from  step  (3)  for  0=0°  and 
for  0=0.  Using  equation  (114),  find  values  of  8  for  0=0°. 

(0)  With  the  aid  of  step  (4)  obtain  approximate  values  of 
Od  associated  with  the  values  of  6  obtained  in  step  (5).  Then, 
with  0=0i,  use  equation  (1)12)  to  obtain  an  approximate 
value  of  0,  where 


_2  in,  F(9) 


0= 


and  the  leading-  and  trading-edge  values  of  0  are  obtained 
from  equation  (B15). 

(7)  Calculate,  from  equations  (Bo)  and  (BlO),  exact  values 
of  ^2 7r  jp.  2ir  jjj,  associated  with  the  initial  values  of  (0,  0)  in 

step  (fl)  where  i=xx-\-x,  and  ?/=?/,  + 

(8)  If,  on  comparison  of  the  coordinates  in  step  (7)  with 
the  coordinates  in  step  (2),  it  is  deemed  necessary  to  approxi¬ 
mate  (0,  0)  more  closely  for  several  of  the  points  (/,  y),  one 

({^ 

procedure  is  the  following:  An  expression  for  /$)  can 

he  found  from  equations  (4)  to  (0)  as 

_ it _ = 

</(0+(0) 

d_  r  sinh  (i p+i0)  sinh  (0-?-f0)  “I 

*2 ir  <)s  ^  |_eosh  7,,+eosh  (0  +  /0A  c osh  70— cosh  (0+f0)J 


■  d. 

- 1  —sin  0 

2  IT 


[cosh  ( 
sinh  7n-f- si 


(0  +  »0) 


sinh  (0+(0)  "^sinh  7,,— sinh  (0+70)  _| 


cosh  (0+/0)  ] 


With  the  notation  of  equations  (Bo)  and  (BlO).  this  expres¬ 
sion  may  he  written 

<sf)  +  . 

iwn»rv+“1 

=  cos  0  sinh  (0  +  701  (  *  e~'*i+  — 

\p  1  Pi  / 

—  i  sin  0  cosh  (1 0  —  70)  (  '  c"l<>;j+  -  c“llt4  ) 

\P3  Pi  ' 

r  .  ,  /  cos  4>\  ,  cos  0>\ 

p  =  COS  0  I  Slim  0  COS  0  1  H - -  —j 

,  ,  •  „  /  sin  0,  ,  sin  02Y1 

+  eosh  0  sin  8  (  - 1 — - — 


where 


V  Pi 


+sin0  [sinh  0  sin  e(^+C-^) 
-cosh  0  cos  0(^+^)] 


and 


,  r  ,  ,  .  .  / cos  0,  ,  cos 
<2= cos  0  cosh  0  sin  0  I - 


_  +■-*) 
\  Pi  Pi  / 


■  1  ,  . /sin  0i  sin  4>.\1 

—sinh  0  cos  0  ^ 

11  . /cos  03  ,  cos  0A 

—  sin  0\  cosh  0  cos  0  ^  — — J 

,  -i  ,  •  a  /sin  03  .  sin  0A1 
—sinh  0  sin  8  (  -  +  --JJ 


The  following  relation  may  then  he  noted: 

YTM2”') 


A0  +  /A01 


/»  +  '"/ 


(BIS) 


Let 


A/M+Y+Y 

a 


where  the  subscripts  0  and  1  refer  to  the  coordinates  given 
in  steps  (2)  and  (7),  respectively.  If  the  values  (0,  0) 
obtained  in  step  ((>)  are  used,  evaluation  of  equation  (B18) 
gives  values  (A0,  A0),  and  (0  +  A0.0  +  A0)  represents  the  next 
approximation  to  the  desired  coordinates.  The  process  in 
steps  (7)  and  (8)  can  he  repeated  if  deemed  necessary. 


APPENDIX  C 

DERIVATION  OF  EQUATION  (12) 


The  transformation  (equation  (0))  from  the  z'-  to  the 
2-plane  may  be  rewritten 

log  ^=f(z) 

=i£  (Cl) 

I  3 

where  the  complex  constants  c„  may  be  defined  as 


Substituting  equation  (('.'!)  in  etpiation  (Cl)  yields 


/(*)=;!  +(«)s 

*  j  ii  i  - 

(C4) 

For 

\R\  .  ... 

i <(  1 ,  thi'  geometric  scries  m  equation 

(C-l)  can  he 

summed  and 

I  n*  Rei: 

(C5) 

Cn=a,+ib„ 

On  the  boundaries,  z'=ae*+,e  and  z-=<ie*n*i‘>\  hence, 


and 


where 


z 

cos  n<t>+j{i  sin  n <*>) 
R=ae*° 


(C2) 


With  0  considered  as  a  function  of  <j>  denoted  by  the 

coeflicients  in  equation  (C2)  are  obtained  us 


^(0)cos n<t>d<t> 

b  1  C2’ 

jg5=-Jo  ^(0)sin n<t>d<t> 


(C3) 


which  can  immediately  be  expressed  as  in  equation  (12). 

For  computational  purposes,  equation  (12)  may  he 
separated  into  real  and  imaginary  parts.  Let  t(z)=p+iq 
and  z = .r  -f-  iy  (where,  for  example,  in  the  smith  approxima- 


wltere 


Then, 

1 

1-2. 

V,  , 

T(0) 

i/O 

77  'h 

1 

2=^. 

A’, 

Jo 

7 7d+ 

X  If 

Ar,  —  jt  cos  sin  0—  1 


X  .  7/ 

S 2~Tl  s*n  ^ ~~Ji cos  & 


D.  l-2(Jcos*+|,i„*)+^ 

and  the  integrations  can  bo  conveniently  performed  by 
Simpson’s  rule. 


APPENDIX  D 

TRANSFORMATION  FROM  r-PLANE  TO  ic-PLANE 


The  linear  fractional  transformation 


w= 


nz+b 

cz+d 


on  which  the  derivation  of  equation  (13)  is  based,  has  the 
following  well-known  properties: 

(1)  When  z  traverses  a  circle  C„  w  traverses  a  circle  C„. 

(2)  Two  points  wx  and  w2  inverse  with  respect  to  a  circle 
C„  correspond  to  two  points  zt  and  z2  inverse  with  respect  to 
the  circle  C\. 

(3)  The  anharmonic  ratio  of  four  points  is  preserved;  that 
is,  if  zx,  z2,  zz,  and  z4  correspond  to  wlt  uy  uy  and  u\. 


(z,  —  zt)(z2—z3)  _  (W|  —  ir,)(w-2—  Hh) 

(z, — z2)  (Zz -  Zi)  (w,  -  uy  (w3  -  wt) 

For  the  desired  correspondence  it  is  known  that  four  points 

£2  _  g2 

w,=6,  w2  =  —b,  and  their  inverse  points  «’3=-p  u’«  =  — £ —  are 

to  correspond  to  s,=/3 1(  z2  =  —  02  and  their  inverse  points' 

R2  —R2 

z3=—>  2,=-^ —  Property  (3)  yields  a  relation  that  mav  be 
01  02 

used  to  solve  for  the  radius  .S'  and  that  can  be  expressed  by 
equations  (15)  and  (10).  When  the  radius  of  the  circle  in 
the  w-plane  has  been  determined,  property  (3)  can  again  be 
used  by  replacing — sav,  wt  by  w  and  z4  by  z.  This  procedure 
will  yield  a  result  that  is  equivalent  to  equation  (13). 
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APPENDIX  E 

OUTLINE  OF  CALCULATION  PROCEDURE 


(1)  List  nirfoil-spction  coordinates  in  percent  chord. 

(2)  Choose  axes  (appendix  B,  paragraph  following  equa¬ 
tion  (Bl7)). 

(3)  List  stagger  angle  0  and  find  y0  and  value  of  c/il  for 
the  skeleton  line  lattice  (table  I,  fig.  2.  and  equation  (8)). 

(4)  Find  (i p,  d)  (appendix  B). 

(5)  Find  t(</>)  (equation  (10))  by  method  given  in  appendix 
of  reference  3. 

((’>)  Blot  i p  against  <t>  where  4>=d^-t.  Find  constant 
(equation  (11))  and  lt=ae*°. 

(7)  Find  complex  constants  0,  and  0,  (equation  (12)  and 
appendix  C). 

(8)  Find  constants  cosh  -y i »  7n  and  /v  =  lq  +  /L  (equa¬ 
tions  (16),  (15),  and  (14)). 

(9)  Find  complex  constant  S  —  'Se1'  (equation  (26))  and 
real  constant  X  (equation  (27)).  Also  obtain  functions 
m(<t>)  and  m(0)  from  equation  (29). 

(10)  Find  a  and,  in  particular,  av,  (equation  (30)). 

(11)  Evaluate  factors  B,  C,  D,  and  E  (equation  (35)). 

(12)  Evaluate  factor  A  in  equation  (35),  first  choosing 
the  angle  of  attack  a  as  indicated  in  the  following  paragraphs: 

The  lift  coefficient  is  as  in  equation  (37) 


r  _ i  _L  _*  _9  2 

L  cjd'lVd  ~  c 


d 

V 


Here  c/d  refers  to  the  value  of  r/d  at  0  percent  chord  minus 
x/d  at  100  percent  chord.  Bv  using  equation  (19)  for  r/2  Vd, 
Or.  may  be  expressed  as 

C’t=//sin  (a-F>i)  (El) 


where 


and 


H-4d- 

c|_\cosh  7i/  \sinh  7i/  J 

.  ( sin  <r„  cosh  7,\ 

an  \cos  o-n  sinh  7,/ 


This  relation  may  be  used  to  find  a  for  auv  desired  value 
of  CL  and  it  is  further  noted  that  a  =  —  17  is  the  angle  of  zero 
lift. 

The  “ideal”  angle  of  attack,  introduced  by  Theodorscn,  is 
defined  for  a  thin  section  as  the  angle  of  attack  for  which  a 
stagnation  point  exists  not  only  at  the  sharp  trailing  edge 
but  also  at  the  sharp  leading  edge.  For  thick  airfoils,  the 
ideal  angle  of  attack  is  defined  in  the  same  manner  (the 
pressure  difference  at  the  hauling  edge  vanishes)  although 
the  point  that  is  considered  the  leading-edge  point  is  not 
precisely  defined.  If  this  point  is  taken  to  be  the  intersection 
of  the  /-axis  with  the  airfoil  leading  edge,  the  ideal  lift  and 
ideal  angle  of  attack  are  found  as  follows:  Let  <q  be  the  value 
of  a  corresponding  to  the  leading-edge  point.  The  quantity 

in  equation  (24)  (or  the  factor  A  in  equation  (35)) 


vanishes  for  The  relation  that  gives  the  value  of 

the  ideal  angle  of  attack  a— a,  is  then 

sin  a _  cosh  y,  sin  a , — sin  <r0 
COS  «  sinh  7i  COS  crL  — -  COS  (To 

and  the  ideal  lift  coefficient,  from  equation  (El),  is 

A  1  1  /  , 

<  lj=~ 4-  jcos  2  (°'i  ao) 

where 

J2=£cosh  7(  cos  |(o-l-t-o'o)l  +  [jsinh  7,  sin  ^(cq  +  ero)  j 

(13)  The  surface  velocity  ratio  v/V  is  now  found  from 
equation  (35)  and  the  (mean)  superstream  pressure  is  found 
from  Bernoulli’s  equation  as 


The  angle  through  which  the  stream  is  turned  may  be  found 
from  equation  (23). 

A  remark  may  bo  inserted  here  regarding  an  inverse 
calculation  procedure.  Instead  of  starting  with  a  given 
lattice,  it  may  be  convenient  to  start  with  given  function 
't'to),  (quantity  i  as  a  function  of  </>)  and  given  parameters 
7o  and  0.  Then  both  the  lattice  arrangement  and  the  flow 
properties  follow  uniquely  and,  in  this  way,  systematic 
families  of  lattices  can  be  studied. 
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Positive  directions  of  axes  and  angles  (forces  and  moments)  are  shown  by  arrows 


Axis 

Force 
(parallel 
to  axis) 
symbol 

Moment  about  axis 

Angle 

Velocities 

Designation 

Sym¬ 

bol 

Designation 

Sym¬ 

bol 

Positive 

direction 

Designa¬ 

tion 

Ssr 

Linear 
(compo¬ 
nent  along 
axis) 

Angular 

X 

X 

L 

-• 

Y - >Z 

Roll 

u 

P 

Lateral  _  .. 

Y 

Y 

Pitching..  .. 

M 

z — >x 

Pitch 

6 

V 

q 

Normal _ _ 

Z 

z 

N 

X - >F 

* 

iff 

r 

Absolute  coefficients  of  moment  Angle  of  set  of  control  surface  (relative  to  neutral 

r  _  L  _  ___M  r  _  N  position),  S.  (Indicate  surface  by  proper  subscript.) 

c,-£bS  °m~qcS  La~c[bS 

(rolbng)  (pitching)  (yawing) 


4. -PROPELLER  SYMBOLS 


D 

V 

Diameter 

Geometric  pitch 

P 

P 

Power,  absolute  coefficient 

p/D 

V' 

Pitch  ratio 

Inflow  velocity 

c. 

Speed-power  coefficient 

v. 

Slipstream  velocity 

V 

Efficiency 

T 

Q 

Thrust,  absolute  coefficient  CT= — 

’  pn2Lr 

Torque,  absolute  coefficient  CQ = p^jy, 

n 

4> 

Revolutions  per  second,  rps 

Effective  helix  angle =tan-1^fj^^ 

5.  NUMERICAL  RELATIONS 

I  hp=76.04  kg-m/s=550  ft-lb/sec  1  lb=0.4536  kg 

1  metric  horsepower=0.9863  hp  1  kg=2.2046  lb 

I  mph=0.4470  mps  I  mi=  1,609.35  m  =  5,280  ft 

1  mps=2.2369  mph  I  m=3.2808  ft 
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